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This project began with the acquisition of short-wavelength, high-dispersion IUE 
spectra of selected late 0- and early B-type stars that are near the main sequence in 
open clusters and associations. The profiles of the resonance lines of N V, Si IV, and 
C IV were studied, and we found that the C IV lines are the most sensitive indicators 
of mass loss (stellar winds) in stars of this type. The mass loss manifests itself as an 
extension of the short-wavelength absorption wing of the doublet, while there is no P 
Cygni-type emission on the long- wavelength side of the line profile. 

We investigated whether the short-wavelength extension could be caused by 
blended lines of other ionic species formed in the photosphere. Although blending 
is present and introduces uncertainty into the estimation of the precise location on 
the main sequence of the onset of the mass-loss signature, it is a crucial issue only 
in a few marginal cases. Mass loss certainly overwhelms blending in its influence on 
the spectrum between spectral types BO and Bl (effective temperatures in the range 
25,000-27,000 K). 

We defined a parameter called P w to describe the degree of asymmetry of the 
C IV resonance-line profile, and we studied the dependence of this parameter on the 
fundamental stellar parameters. For this purpose, we derived new estimates of the 
stellar T e ff and log g from a non-LTE, line-blanketed model-atmosphere analysis of 
these stars (Grigsby, Morrison, and Anderson 1992). In order to estimate the stellar 
luminosities, we performed an exhaustive search of the literature for the most reliable 
available estimates of the distances of the clusters and associations to which the 
program stars belong. 

Thus, we were able to study the dependence of P w on stellar temperature and 
luminosity. P w shows a sharp increase at a luminosity of about 10 44 L©, while the 
total equivalent width of the C IV doublet increases more gradually, since these 
lines are present in the photospheres of stars that Eire below the wind theshold. As a 
function of T e ff, the total equivalent width increases gradually across the threshold, 
and, below the threshold, it is well matched by the plane-parallel, line-blanketed, 
non-LTE model atmospheres of Grigsby, Morrison, and Anderson (1992). This good 
agreement indicates that, in this temperature and luminosity regime, additional 
ionization sources such as X-rays are not needed to account for the strength of the 
C IV doublet. 

The temperature and luminosity threshold at which we first find stellar winds 
corresponds well to an effect in the population of the C IV ground state that 
occurs in our model atmospheres at T e ff = 28,000 K. Below this temperature, the 
population of this state is smaller at the top of the photosphere than it is in deeper 
layers. At this and higher temperatures, however, the ground-state population is 
increasingly enhanced relative to that in deeper layers. This population enhancement 
is a non-LTE effect. Although we have not investigated the consequences of this 
effect in detail, we speculate that perhaps radiation pressure on this enhanced 
population of C +++ ions could initiate a radiation-driven wind at this temperature. 

This research employed the SIMBAD database, which is operated by the Centre 
de Donnees Astronomiques de Strasbourg (France). This research also comprised 
part of the Ph. D. research of James A. Grigsby (U. of Toledo). 
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ABSTRACT 

The profiles of ultraviolet resonance lines of C rv were surveyed in a sample of 29 cluster and assocation 
members in the spectral type range 09-B2 III-V, together with a few field stars of interest. The temperatures 
and gravities of the stars were taken from the model atmosphere analysis by Grigsby, Morrison, & Anderson 
(1992), and the luminosities were estimated on the basis of cluster and association distances from the recent 
literature. A parameter P w , was defined in order to describe the degree of asymmetry of the C iv profile. This 
parameter, together with total C iv equivalent width, was found to be well correlated with stellar luminosity 
and temperature. A few anomalous stars were noted: r Sco, HD 66665, HD 13621, and the ON stars 
HD12323 and HD 201345. The results suggest a sudden onset of observable mass loss at T cf( = 27,500 ± 500 
K, log (L/L q ) = 4.4 ± 0.1, in agreement with the previous study by Prinja (1989). At T rff = 28,000 K and log 
0 — 4, our non-LTE model atmospheres show an enhancement in the ground-state population of C +3 in their 
topmost layer, which could be responsible for initiation of the winds via radiation pressure on the C +3 ions, 
or for the onset of visibility of C* 3 ions in the wind because of an increase in the optical depth in the C iv 
lines in the outermost layers. 

Subject headings: stars: early-type — stars: mass loss — ultraviolet: stars 


1. INTRODUCTION 

This paper is a survey of ultraviolet resonance line profiles in 
a sample of main-sequence stars in the spectral range from 09 
to B2, which straddles the theoretical limit of stability of these 
stellar atmospheres against radiation pressure (Abbott 1979). 
Our purpose in this survey is to determine more precisely the 
location on the main sequence of this limit and the dependence 
of the morphology of the C iv resonance line profile on funda- 
mental stellar parameters above this limit. 

Previously, most such studies have concentrated on more 
luminous stars with stronger winds (e.g., Pauldrach et al. 1990). 
Studies of O stars (Garmany et al. 1981; Garmany & Conti 
1984) found that mass-loss rates in these stars depended upon 
luminosity roughly in the manner predicted by the radiation- 
driven wind theory (Castor, Abbott, & Klein 1975; Abbott 
1982). but with much unaccounted-for scatter. Investigations 
of the C iv /U1548, 1550 resonance lines in B stars (Barker 
1984) showed a strong correlation between the C rv line 
strength and spectral type, with a large scatter at each spectral 
type; Walbom & Panek (.1984) found a similar trend — without 
the scatter — on the basis of a smaller sample of stars. Finally, 
some main-sequence stars have been identified by Walbom, 
Nichols-Bohlin, & Panek (1985) as having enhanced winds for 
their spectral types. 

More recently, the behavior of wind signatures in this range 
of spectral types and luminosity classes has been studied by 
Grady, Bjorkman, & Snow (1987), Howarth & Prinja (1989), 
and Prinja (1989). Our study differs from theirs in at least two 
ways. First, all of them studied the wind diagnostics as a func- 
tion of spectral type, or else in terms of a temperature scale that 

* Guest Observer with the WE satellite, which is operated jointly by 
NASA. ESA, and the SERC. 


is calibrated in tenns of spectral type, leaving open the possi- 
bility that some of the observational scatter that they found is 
due to errors or to the binning effects in spectral classification. 
In contrast, we have studied the mass-loss signatures of our 
sample of stars as a function of effective temperature and 
gravity. In addition, our non-LTE, line-blanketed model atmo- 
spheres provide more insight into the physical conditions in 
the photosphere underlying the wind. Second, they studied a 
larger range of spectral type and/or luminosity class. Here we 
have used a sample of main-sequence stars that covers the 
region in the H-R diagram around the static limit at least as 
densely as the samples used in those investigations. 

This is the second of two papers on late O and early B-type 
stars near the main sequence. In the first paper (Grigsby, Mor- 
rison, & Anderson 1992, hereafter Paper I) we reported on the 
analysis of a sample of stars with a line-blanketed, non-LTE 
model atmosphere code (PAM). In that analysis we derived 
effective temperatures, surface gravities, and projected rota- 
tional velocities from hydrogen and neutral helium lines, and 
we studied the behavior of the equivalent widths of N n, N m, 
C n and C m lines as a function of effective temperature in a 
search for abundance differences from star to star. 

In the present paper, we examine the morphology of wind- 
sensitive ultraviolet resonance lines — N v 221239, 1243, Si rv 
221394, 1403, and C rv 221548, 1550 — and we conclude that 
the best probe of onset of the stellar winds is the C iv doublet 
We define the “wind parameter” P w , which measures the 
asymmetry of this doublet and we discuss its dependence and 
that of C tv total equivalent width on stellar temperature, 
gravity, and luminosity. 

Of course, it is conceivable that some stars have winds that 
are too transparent to be observed even in this sensitive spec- 
tral diagnostic, and it is widely recognized that some of the 
material in the wind is likely unobservable in any case because 
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it is ionized by shocks to states higher than C* 3 . A good 
discussion of the relevant literature on this point is given by 
Prinja, Barlow, & Howarth (1990). In hot stars. X-ray emission 
is thought to be produced in shocks that are created by insta- 
bilities in radiation-driven winds (e.g^ MacFarlane & Cassin- 
elli 1989), and this emission may be an even more sensitive 
diagnostic of mass loss. Recently, X-ray surveys of main- 
sequence stars (Cassinelli et al. 1994; Meurs et ai. 1992) have 
shown that X-ray emission persists farther down the main 
sequence than does asymmetric C rv absorption. The available 
evidence is not, however, sufficient to indicate precisely how 
much farther in terms of effective temperature. One goal of this 
investigation is to locate more precisely the threshold in the 
theoretical H-R diagram where stellar winds first appear in the 
C iv doublet and to identify a mechanism that could cause 
either the initiation of the wind itself or the development of 
sufficient optical depth in the wind to cause asymmetric C rv 
absorption. 

In the following, § 2 describes the observational material, § 3 
discusses the stars’ luminosities and, thence, locations in the 
H-R diagram, and § 4 discusses the S rv and N v lines. In § 5, 
we describe our efforts to disentangle the effects of blends on 
the C iv line profiles, we define the wind parameter, and we 
discuss its dependence and that of C iv equivalent width on 
stellar temperature and luminosity. Section 6 discusses C rv 
line formation in our static, non-LTE, line-blanketed model 
atmospheres, and § 7 is a discussion of the initiation of winds 
in light of our results and those of recent X-ray observations of 
hot stars. Section 8 presents our conclusions. 


2. OBSERVATIONAL MATERIAL 

The program stars are described more fully in Paper I, along 
with their selection criteria. Ouster membership and other 
information may be found in Table 1 of Paper I. The star HD 
143275 (<5 Sco) was included, even though it has not been care- 
fully checked for photometric constancy, because its angular 
diameter is known; the star HD 66665, though not a cluster or 
association member, was included because Massa (1985) has 
identified that star as one with unusually strong C tv resonance 
lines. Some OBCN stars have also been included in violation of 
the same requirements because of their abundance anomalies 
(Walborn 1976). 

Table 1 of the present paper lists the stars from Paper I for 
which WE observations exist The data were processed with 
the standard procedures (Tumrose & Thompson 1984). 
Observations made by us were processed immediately; the 
remaining images, which were taken from the IUE archives, 
were processed in early 1987 with the then-current software. 
All the spectra are well exposed in the region of the C iv 
doublet LAI 548, 1550, but the spectra of the more highly 
reddened stars tend to be less well exposed at shorter wave- 
lengths. Thus, the data quality is typical of WE spectra, with 
signal-to-noise ratio similar to that estimated for individual 
images by Leckrone & Adelman (1989). The unsmoothed 
spectra have the nominal spectral resolution (Grady & Imhoff 
1985). The data were continuum normalized using the routine 
IUESPEQ developed by the WE Regional Data Analysis 
Facility (RDAF) staff at Goddard Space Flight Center (IUE 


TABLE I 


Luminosities of the Program Stars 


Star 

(HD/HDE) 

Spectral Type 

SWP 

Association 

Subgroup/Cluster 

Reference 

5 log r — 5 

A(5 log r— 5) 

log (LJL 0 ) 

A log (LJLq) 

12323 

ON9.5 V 

9652 

Per OBI 

1 

11.8 

0-5 

4J90 

020 

13621 

B0.5 IV 

30483 

Per OBI 

1 

11.8 

0.5 

4.896 

020 

21856 

B1 V 

20970 

Per OB2 

2 

7.5 

0.2 

3.760 

013 

24131 ....... 

B1 V 

20845 

Per OB2 

2 

7 3 

0.2 

3.870 

013 

34816 

B0.5 V 

26726 

On OBI 

3 

7.9 

0.6 

4.450 

026 

35299 

B15 V 

18005 

Ori OBla 

3 

7.9 

as 

3.780 

022 

35653 

B0.5 V 

30479 

Aur OBI 

4 

10.5 

0.4 

4-578 

0.14 

36430 

B2 V 

29158 

Ori OBlc 

3 

7.9 

a6 

3.430 

026 

36822 

B0 IV-V 

29099 

Cr 69 . 

5 

8.0 

0.2 

4.540 

0.16 

36824 

B2 V 

29165 

Ori OBla 

3 

7.9 

0.5 

3-260 

0.22 

37209 

B1 V 

29160 

Ori OBlc 

3 

7.9 

0.6 

3.786 

026 

37366 

09.5 V 

30480 

Aur OBI 

4 

10.5 

0.4 

4.650 

0.14 

120307 

B2 IV 

36482 

Upper Cen 

6 

5.7 

03 

3.780 

019 

132058 

B2 III 

17458 

Upper Cen 

6 

5.7 

03 

4.030 

019 

142669 

B2 IV-V 

22091 

Upper Sco 

6 

6.0 

OJ 

3.620 

016 

142375 

B0.3 IV 

17395 

Upper Sco 

6 

6.0 

0.3 

4.590 

015 

144470 

B1 V - 

13609 

Upper Sco 

6 

60 

OJ 

3560 

015 

148703 

B2 III 

29161 

Upper Sco 

6 

60 

OJ 

3.610 

019 

149438 

BO V 

13793 

Upper Sco 

6 

60 

0-3 

4J90 

016 

151985 

B2 IV 

5879 

Upper Sco 

6 

60 

OJ 

3.750 

016 

152246 

09 III 

16214 

Sco OBI 

7 

IIJ 

OJ 

5.410 

014 

163892 

09 IV 

29157 

Sgr OBI 

8 

. IIJ 

0.1 

5.160 

013 

207538 

BO V 

8992 

Cep OB2 

1 

9.9 

03 

4.850 

017 

208440 

B1 V 

30482 

NGC 7160 

1 

9.9 

OJ 

4.070 

019 

214680 

09 V 

6784 

Uc OBI 


8.5 

0.5 

4.640 

021 

217101 

B2 IV-V 

18309 

Lac OBI 

1 

65 

OJ 

3.830 

021 

218195 

09 III 

26975 

Cep OBI 


11.8 

OJ 

5J70 

025 

218342 

BO IV 

29164 

Cep OB3 


9J 

03 

4.680 

016 

251847 

B1 IV 

30481 

Gan OBI 

9 

11.8 

0.1 

4.570 

015 


References.— (I) Garmany &. Stenccl 1992. (2) Strom. Strom. & Carrasco 1974; IC 348. (3) Anthony-Twarong 1982. (4) Kimeswenger & Weinberger 
1989; NGC 1960. (5) Murdin & Pension 1977. (6) Jones 1971, averages quoted in de Gcus et aL 1989. (7) Perry et aL 1991. (8) McCall, Richer, & 
Visvanathan 1990; NGC 6523. (9) Chavama-K, de Lara, & Hasse 1987; S252/NGC 2175. 
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RDAF User’s Tutorial Manual 1990). This software performs 
ripple corrections (Barker 1984) automatically and allows the 
user to view and normalize large segments of spectra. The 
spectra were normalized by fitting straight line segments to 
high points across regions several hundred angstroms wide 
and were smoothed with a five-point boxcar smoothing pro- 
cedure. The continua produced in this manner by several inde- 
pendent but experienced users of the software are generally 
astronomer independent to within a few percent and agree well 
with the predictions (Grigsby 1991 ; see also Figures 5a and 5b 
of this paper) by the non-LTE, line-blanketed model atmo- 
sphere code described in Paper I. 

3. THE PROGRAM STARS IN THR H-R DIAGRAM 

Table 1 lists the luminosities of the duster and assodation 
members and the estimated uncertainties in the luminosity. 
The B — V colors were taken from Humphreys (1979), Bcrtiau 
(1958), and the photometric references found in Table 1 of 
Paper I. The intrinsic colors used were those of Johnson (1966), 
and A, was computed with R v [the ratio of the total absorption 
in the V band to the color excess E(B — V)] equal to 3.1, follow- 
ing the standard reddening law. Bolometric corrections were 
taken from Buser & Kurucz (1978), who obtained them by 
fitting model (LTE) stellar atmospheres to observed fluxes in 
main-sequence B stars. We have adopted the additive correc- 
tion of 0.1 to all the bolometric corrections that those authors 
recommend. 

We searched the literature for the best available information 
on the distance moduli of OB associations, and the results are 
incorporated into Table 1. For the Sco-Cen assodation, the 
subgroup assignments are taken from de Geus, de Zeeuw, & 
Lub (1989), and for Ori OBI from Warren & Hesser (1978). 
For distance moduli, we adopted values from the recent review 
by Garmany &. Stencel (1992), the reanalysis of Ori OBI by 
Anthony-Twarog (1982), and the reanalysis of Sco OBI by 
Perry, Hill, & Christodoulou (1991). For assodations not 
included in those references, we could find no recent estimate 
of the distance to the assodation as a whole. Therefore, we 
used the best available estimate of the distance modulus of the 
open cluster belonging to each assodation. In those cases, the 
cluster is named in the list of references to Table I. For some 
clusters, such as NGC 2175, the distance moduli in the liter- 
ature differ by as much as 1 mag. In these cases, we chose the 
estimate that we judged to be based on the soundest method- 
ology and/or the largest number of cluster members. 

The uncertainties in luminosity were computed from the 
sum, in quadrature, of the uncertainties of the quantities that 
constitute the bolometric magnitude: the bolometric correc- 
tion, m„, A v , and the distance modulus. The estimated errors in 
these quantities, except for the distance moduli, are upper 
limits. For the errors in the distance moduli, either we used the 
uncertainty quoted by the original investigator, or, if no error 
was quoted, we estimated the uncertainty from the scatter of 
the main-sequence stars in the cluster H-R diagram. In addi- 
tion, we estimated the uncertainty in the distance modulus that 
is due to the radial extent of the cluster or association by 
assuming that the object in question is spherical with a radius 
computed from the angular diameter and distance. These 
errors in the distance modulus are lower limits in that they 
include only the front-to-back uncertainty in a given cluster 
(association) due to its finite size, but not any systematic errors 
that may be present. In modem estimates of distance moduli, 
the formal error in the mean distance modulus is often smaller 


than the standard error of the individual stars' moduli about 
the mean, which may indeed, therefore, be dominated by front- 
to-back scatter in the stellar distances. In our analysis, we took 
the uncertainty in the distance modulus of a star to be the error 
quoted for the association or the front-to-back uncertainty, 
whichever is larger. The uncertainty in A v results from errors in 
R„ and the color excess E(B — V). We took the uncertainty in 
R , to be 0.2 and the relative uncertainty in E(B — V) to be 15%, 
an extreme upper limit. Finally, we took the uncertainty in the 
apparent visual magnitude to be 0.05 and that of the bolo- 
metric corrections to be 0.2; both are conservative upper 
limits. Even with these upper limits, the uncertainties in the 
luminosities in Table 1 are generally dominated by the uncer- 
tainties in the distance moduli, except for stars that are heavily 
reddened (A, > 0.7), in which case the error in A„ is compara- 
ble to or larger than that in the distance modulus. 

For binaries in which the V magnitudes of the secondary 
components are known and in which the separations are less 
than 10", we corrected the combined V magnitude to that of 
the primary alone. For larger separations, we assumed that the 
published V magnitude refers to the primary alone, and we 
made no correction. Thus, small corrections were needed for 
HD 13621 and HD 37209, for which we obtained the magni- 
tudes and separations from the HIPPARCOS Input Catalog 
(1992). For HD 143275 (6 Sco), we corrected for the companion 
at A m - 1.5 observed by Bedding (1993). 

The H-R diagram is shown in Figure 1. Included are the 
theoretical zero-age main sequence (ZAMS) ( solid line) and the 
5 x 10 6 year isochrone ( dash-dotted line ) from Maeder & 
Meynet (1988). Since this age is typical of our program stars, 
their location just to the right of this isochrone is well 
explained by stellar evolution. Also included in the figure is the 
main-sequence boundary at the low-temperature extreme 
given in Maeder & Meynet (1988). Although our sample of 
stars was not chosen with the goal of testing the boundaries of 
the main-sequence band in mind, the sample is indeed confined 
to the theoretical main sequence. 

An examination of the surface gravities derived in Paper I 
reveals some minor discrepancies. For example, at an effective 
temperature of 28,000 K, the lower gravity star HD 36822 
apparently has a lower luminosity than HD 13621; at T c(( = 
22,000 K, HD 217101 (log g = 4.2) apparently has a higher 
luminosity than several stars with significantly lower gravity. 
However, these discrepancies are all easily accounted for by the 
uncertainties in the luminosities, and none is significant. 

4. EVIDENCE FOR WINDS IN THE N v AND Si IV LINES 

Walbom et al. (1985) discuss this topic and display data for 
some of the stars in the present sample. They find evidence for 
winds in the N v <111239, 1243 resonance lines in the so-called 
nitrogen-enhanced stars HD 12323 and HD 201345, and in 
t Sco, which they label as peculiar. These authors also suggest 
that the Si rv resonance lines in r Sco display weak evidence for 
a wind through their slight asymmetry. 

Figures 2a and 2b illustrate the evidence for winds in the N v 
resonance doublet in selected stars in the present sample. The 
stars HD 218195 and HD 163892 may also show evidence for a 
wind in N v, but both spectra are extremely noisy due to 
having been exposed incorrectly and have been omitted. In this 
group, the lines in HD 12323 and HD 201345 are not asym- 
metric, but are strong and shortward shifted. The lines in HD 
214680 are weak, and any asymmetry seen may be the result of 
blends — the same is true of HD 37366 and HD 218342 (not 
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Effective Temperature (K) 

Fig. 1. — H-R diagram of the program stars. The lines are from Maeder & Meynet (1988). The solid line is the theoretical ZAMS. The dash-dotted line is the 
5 x I0 6 yr isochrone, and the dashed line is the low-temperature boundary of the main sequence. Open drdes refer to stars with log g S 3.8; solid dots, log g =■ 4.0; 
open squares, log g > 4.0. 


shown). In HD 152246 we see evidence of a broad, shortward- 
shifted absorption trough, with evidence of P Cygni-type emis- 
sion, though the noisiness of the spectrum makes the latter 
conclusion uncertain. In the present sample, the most striking 
N V lines remain those of t Sco, but HD 66665 emerges as a 
dose second. Its N v lines are stronger and more asymmetric 
than those of HD 12323, HD 201345, and HD 214680, even 
though all of the latter stars have higher effective temperatures 
according to Paper I. In general, HD 66665 displays unusually 
strong resonance lines, as shown below and in Grigsby (1992). 
No longer should t Sco be viewed as unique for its spectral 
type (Walbom et al. 1985). 

Figures 3a and 3 b show stars that exhibit asymmetries in 
their Si rv resonance lines. The most pronounced asymmetries 
occur in HD 66665, HD 152246, and r Sco. The asymmetry 
that appears in HD 207538 and HD 218342 may be due to a 
blend in the wing of the shortward component. The longward 
component is contaminated by a reseau. Once again, HD 
66665 stands out: among the program stars, it has by far the 
strongest and most asymmetrical Si rv doublet. 

Clearly, evidence for winds does occur in the N v and Si tv 
doublets in some of the program stars. Nevertheless, these lines 
are difficult to use in the study of wind development because of 
the weakness of the lines themselves, the weakness of the wind 
signature, or both. Part of this problem is the nature of the 
lines and the location (in temperature) of the H-R diagram 
where the winds switch on: the ionization balance is only 
beginning to switch to N v in the hottest of the program stars. 
Conversely, the ionization balance is shifting away from S iv in 
those stars. The weak dependence of the equivalent width of 
the Si iv doublet on effective temperature (Grigsby 1992) indi- 
cates that this doublet may be an unwise choice for the study of 
weak winds, at least as a means to test the winds’ temperature 
sensitivity. However, it is the weakness of the wind signatures 
in these lines that makes them unsuitable. For the study of 
weak winds, the C tv resonance lines are clearly the best choice. 


5. C iv MORPHOLOGY IN THE PROGRAM STARS 

5.1. Contribution of Fe Blends in the Civ Region 

H'ubeny, Stefle, & Harmanec (1985a, b) and Barbier & 
Swings (1979) have suggested that, in marginal cases, Fe blends 
in the region of the C rv resonance lines may mimic mass loss. 
Because we are interested in determining the precise tem- 
perature on the main sequence where stellar winds first 
produce asymmetries in the C rv profiles, we attempted to take 
blends quantitatively into account As a first step toward this 
goal, we calculated synthetic spectra with PAM, which does 
not include winds, hoping that a close match between the syn- 
thetic and the observed spectrum would allow us to decide 
which features are blends and which are due to wind-enhanced 
C iv. Figure 4 shows the synthetic spectrum produced by PAM 
superimposed upon the observed spectrum of HD 37209 
( T t ff = 25,000 K, log g = 4.0; Paper I), a spectrum that appears 
to be affected by photospheric blends only, not by a wind. The 
synthetic spectrum was calculated by the methods described in 
Paper I, except that all available opacity sources in the 1525— 
1560 A interval (from Kurucz 1987) were included in the solu- 
tion of the equation of radiative transfer (Grigsby 1991). The 
available opacity sources in the case of iron include only lines 
of Fe n and Fe ni. The resolutions of the synthetic and 
observed spectra are identical. The synthetic spectrum has 
been instrumentally broadened with a Gaussian profile with an 
FWHM of 0.13 A to simulate 1UE instrumental broadening 
(Grady & Imhoff 1985). It has also been rotationally 
broadened with a Doppler profile corresponding to a projected 
rotational velocity of 50 km s " 1 (Paper I). 

It is clear that the model matches many of the features 
(primarily Fe ni) in the region only marginally at best. Many of 
the lines in the Kurucz data possess incorrect wavelengths, and 
many are in error in oscillator strength as well. The “lines" 
(actually, these features are composed of hundreds of individ- 
ual lines) that are reproduced best are those between 1530 and 
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Fig. Z — Evidence for winds in the N v doublet for selected stars. In HD 132246 (09 III), the absorption appears shallow and broad extends off the scale of the 
graph. There is the possibility of emission at 1243.3 A. Note the similarity of the profiles in t Sco and HD 66663 (B0.5 IV). In all cases the HD number of the star 
appears directly above its spectrum, and the spectra are shifted vertically by an additive constant Reseaux and other band data Bagged by the 1VE software are 
indicated by X*s. The vertical lines indicate the rest wavelengths of the N v doublet The apparent redshift of the 212428 line is unexplained. The star's radial velocity 
of 20 km s' 1 does not account for it 


1535 A. The C iv doublet is matched reasonably well, and the 
blends atop the carbon lines also appear in the model but are 
too strong. At 1546.5 A, a line clearly appears in the observ- 
ation which is not predicted by the model. A search of the 
Abbott (1985a) line lists reveals aCm line at this wavelength; 
the same line does not appear in the Kurucz data. It also 
appears in other stars in the sample, such as HD 24131 (T rff = 
24,500 K, log 0 = 4.0). 

If this feature is a line of C in, then it might be expected to 


increase in strength to a maximum at an effective temperature 
of — 28,000 K (see Fig. 1 1 of Paper I) and so may be expected 
to contribute to the asymmetry that we might otherwise have 
attributed to the beginnings of a wind in stars such as HDE 
251847 (T eff = 28,000 K, log g = 4.0). Figure 5b shows what 
might be interpreted as a strengthened line of C in. However, if 
that is the case then the line is redshifted by 0.5 A and asym- 
metric in the redward direction. It is more likely that the line 
seen in HD 37209 (7^ f = 25,000 K, log g = 4.0) is also seen in 
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Fig. 3. — Evidence for winds in the Si rv resonance doublet in selected stars. Note the great strength of these lines in HD 66665 (B0.5 IV). Vertical lines again 
indicate the position of the rest Wavelengths of the doublet 


HDE 251847 at the same wavelength (1545.6 A). This weaker 
blend does contribute somewhat to the asymmetric appear- 
ance of the line, but the main contributor to the asymmetry 
seems more likely to be a blueshifted C iv line. At any rate, the 
blend cannot account for the asymmetry in the longward com- 
ponent of the C tv doublet. 

As a final test of the possibility that iron lines may help 
mimic mass loss, we examine the temperature dependence of 
the blends (Figs. 5a and 56). Published spectra of stars in this 
temperature range (Rountree & Sonnebom 1991; Walborn et 
aL 1985) show 'that the strengths of the blending lines in this 


spectral range vary little with temperature. On the other hand, 
in order to mimic the broadening effect of a stellar wind, a 
blending line must be strongly temperature dependent in the 
temperature range around 28,000 K. Therefore, it is unlikely 
that these numerous Fe lines contribute to this effect To show 
this point more clearly, in Figures 5a and 5b we compare, the 
C rv region in HD 66665 (T^, = 28,000 K, log g = 4.0) and 
HDE 251847 with that in HD 35299 (T cff = 24,000 K, log 
g = 4.0), all stars with relatively low (^20 km s -1 ) projected 
rotational velocities. Since most of the blends here are due to 
Fe in (with the exception discussed above), one would expect 
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FlO. 4. — Comparison of the observed spectrum of HD 37209 ( solid line ) with the synthetic iron lines and C tv doublet from the model (T„ r => 25,000 1C, log 
g « 4.0) presented by Grigsby et al. (1992) ( dashed line). 


that a strong temperature dependence in the lines immediately 
bordering the C iv doublet would be mimicked by similar 
behavior in nearby Fe in lines if the shape of the C rv doublet 
were due to iron blending. Clearly, the lines in the 1540-1545 A 
region are not affected much by the increase in temperature, 
while the C iv doublet is clearly increasing in strength. Indeed, 
the Fe m lines in HDE 251847 appear to get stronger with the 
change in temperature than those in HD 66665, while the C rv 
doublet is clearly weaker in the former star. Fe in lines in the 
1520-1540 A region show behavior similar to those in the 
figures. 

It is clear from the foregoing discussion that the effect seen in 
the C rv lines is not likely to be due to blending due to iron or 
carbon. The most straightforward explanation is that most of 
the asymmetry at 28,000 K is due to a weak wind. 

5.2. Wind Parameter 

The C iv 221548, 1550 resonance doublet is the most sensi- 
tive measure of mass loss in the program stars. Since these stars 
have weak winds, the classical P Cygni profile is not developed, 
and the wind is revealed primarily through the blueward asym- 
metry of the C iv resonance lines. Although actual measure- 
ments of mass-loss rates would be extremely useful in 
understanding the mechanisms which initiate winds, we have 
concluded that we do not have a reliable means of measuring 
the rates. The approach- outlined by Olson (1982) does not 
work for our stars because evidence for winds typically occurs 
only in C tv. Alternatively, we could estimate the mass-loss rate 
from the C rv doublet if the wind’s carbon ionization balance 
were known. With respect to this approach, Lamers & Groene- 
wegen (1990) give empirically derived ionization ratios for C in 
and C tv at half the wind terminal velocity in giants and super- 
giants, but not for main-sequence stars in our temperature 
regime. Therefore, in this paper we do not attempt to derive 
mass-loss rates, but rather characterize the wind-sensitive 
feature with a model-independent parameter, which we then 
examine as a function of the known physical properties of the 
program stars. 


Figures 6a- 6c show the C rv 221548, 1550 resonance doublet 
in three of the program stars. These spectra illustrate the pro- 
gression of the strength and asymmetry in what seems to be a 
fairly smooth sequence. In order to investigate C iv morphol- 
ogy in the program stars we defined a parameter P w that we 
term the “wind parameter,” as follows: 

K = (A2> - A2,)d , (1) 

where d is the maximum depth of the profile relative to the 
continuum and where A2 6 and Al, defined as follows: 

A2> = 2> — 2° and A2 r = 2? - 2, . 

Here, 2* is the wavelength at which the shortward component 
of the doublet intersects the continuum, and 2° is the rest 
wavelength of the shortward component Similarly, 2, is the 
wavelength at which the longward component intersects the 
continuum and 2? is that component’s rest wavelength. Note 
that it is possible for A2 r to be negative if the Doppler shift of 
the doublet is sufficiently large. As an example, Figure 7 illus- 
trates the quantities A2 4 and A2 r in the star HD 207538. These 
definitions ensure that P w is a measure of the asymmetry and 
the depth of the profile. 

Since blends with weak photospheric lines blur the point at 
which a typical wind-sensitive profile intersects the continuum, 
it was necessary to estimate the intersection points. We accom- 
plished this by drawing a profile by hand and extrapolating the 
observed profile to the continuum where necessary. In stars 
with low (^50 km s" 1 ) projected rotational velocities, the 
extrapolation was minimal. As an example. Figure 7 shows the 
hand-drawn profile for the star HD 207538. The circles show 
the points used to define the profile and the point where the 
hand-drawn profile intersects the continuum. One may argue 
that the hand-drawn profile may bend over and intersect the 
continuum at a slightly shorter wavelength; we have chosen 
the profile that is shown because P Cygni profiles generally 
exhibit sharp shortward edges (Prinja et al. 1990). Clearly, 
some ambiguity in the location of the intersection of the profile 
with the continuum remains, and we have attempted to take 
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Fig. 54 

Fig. 5.— Comparison of the C rv region in HD J 5299 (T^, =» 24,000 K, log g - 4.0) with HD 34816 and HD 66665 (both with T M - 28,000 K. log g - 4.0). The 
iron lines show little dependence on temperature, but the C iv lines respond strongly to the temperature increase. 

that ambiguity into account in the estimates of the uncer- absorption: included are contributions due to shortward- 
tainties of the intersection points. shifted discrete absorption components (DA Cs) caused by C rv. 

These hand-drawn points are input interactively (using soft- but contributions from weak photospheric iron lines that 

ware written by J. A. G.) and displayed on a Sun workstation blanket this part of the spectrum are excluded. Like the 

running PGPLOT (Pearson 1990). The dashed line is a result extrapolation of profiles described above, this is a procedure 

of interpolation between the two points, with typically 10 that involves some ambiguity. In the slowly rotating stars it is 

evenly spaced interpolations per interval. The software auto- fairly straightforward to distinguish the carbon absorption 

matically measures the parameters in equation (1) from the from the iron blends (see above), because of the velocity separa- 

input profile and computes the wind parameter. Also com- tion of the discrete C iv components. These stars may be used 

puted is the equivalent width of the input profile. The hand- to make the same identification in the more rapid rotators, 

drawn profile is intended to be a measure of total carbon Since we considered the DACs to be part of the total carbon 



xn U 9*001*8 any aAipjay 


ORIGINAL P^pE IS 

OF POOR QUALITY 








Wavelength (X) 

Fig. 7. — Example of a hand-drawn profile for HD 207538 used to compute the wind profile. Circles denote the location of points chosen by the user, and the 
profile used to determine the equivalent width is delineated by the dashed line. The dashed vertical lines indicate the rest wavelengths of the C tv doublet. The 
horizontal line shows the location of the continuum. The solid vertical line indicates the depth parameter d, and its wavelength represents the wavelength of 
maximum depth. 


absorption, in stars in which the DAC is the deepest part of the 
profile (e.g^ HD 35653) the DAC depth is taken as the depth 
parameter d. 

Nevertheless, it is difficult to separate the effects of Fe and 
C rv absorption in many cases. In the example shown in Figure 
7, the feature that is identified in the longward part of the 
hand-drawn profile as due to carbon may indeed be an iron 
blend. Such misidentifications have no effect on the value of the 
wind parameter, but will affect the equivalent width measure- 
ment. We attempt to account for possible misidentification of 
wind features with a large estimate of the equivalent width 
uncertainty. 

Table 2 shows P w and total C rv equivalent width in the 
program stars, along with the measured quantities A2 S , A2,, 
and d and their estimated uncertainties. The error in P w is due 
to uncertainties in continuum placement and in the measured 
quantities from which it is calculated; since we are fairly con- 
fident of continuum placement, we take the. latter to be the 
dominant contributor, and the uncertainty in P w is the sum of 
the errors, in quadrature, of the errors in A/*, A2,, and d. The 
uncertainty in equivalent width is more difficult to quantify 
and is prone to systematic error as discussed above if photo- 
spheric blends are misidentified as C rv wind features. The 
uncertainties are estimates based on the scatter resulting from 
numerous measurements of several spectra. 

5.3. T emperature and Luminosity Dependence of P w and C rv 
Equivalent Width 

Figure 8 displays the temperature dependence of P w . The 
stars with log g = 4 show a clear monotonic increase in P w 
with effective temperature. If the asymmetry of the C rv lines is 
an accurate indication of stellar winds, then mass loss begins at 
effective temperatures of 26,000-27,000 K for a surface gravity 
of 4 dex, where values of P w become significantly greater than 
zero. Certainly at 28,000 K, the doublet is unquestionably 
asymmetric, but at 25,000 K and below it appears photo- 


spheric (where it appears at all), possessing neither an asym- 
metry nor a radial velocity displacement relative to other 
photospheric lines. DACs do not play a significant role in the 
identification of the critical temperature and gravity where 
asymmetries first appear, simply because none of the stars near 
that region of the H-R diagram display them. The stars that do 
exhibit DACs generally are hotter than 28,000 K or display log 
g < 4.0 (e g., HD 35653, HD 152246, HD 163892, and HD 
218195). As a consequence of the symmetry of the doublet in 
the stars with effective temperatures at or below 25,000 K, P w is 
consistent with zero in those stars and they display no evidence 
of mass loss. 

Two stars with a similar surface gravity that do not agree 
with this picture are the OBCN stars HD" 12323 and HD 
201345 ( T t „ = 29,000 K). Both display weak carbon lines for 
their effective temperature. A possible explanation is that these 
stars are carbon underabundant (Walborn 1976), consistent 
with the results of Paper I, where we found the optical carbon 
lines in these stars to be weak with respect to those in stars 
with similar effective temperatures. HD 12323 shows very 
weak C iv with no trace of asymmetry, while HD 201345 
exhibits a stronger line and only the beginnings of asymmetry. 
The fact that these stars do not display evidence for a wind in 
the carbon lines does not necessarily indicate that they do not 
lose mass. As shown in Figure 2, the N v resonance lines 
(221239, 1243) in HD 12323 appear blueshifted (Walborn et al. 
1985), and in HD 201345 the same lines are blueshifted and 
possibly asymmetric. 

The other stars that do not lie near the log g = 4 locus 
possess significantly lower gravities. This effect is particularly 
apparent in the star HD 35653. At an effective temperature of 
25,000 K, it displays a strong asymmetry in C tv, with obvious 
narrow, shortward-shifted DACs. Stronger wind signatures in 
lower gravity stars are expected because of the lower escape 
velocity and the shift of the carbon ionization balance toward 
higher stages of ionization in the photospheres of such stars 
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Star 

(HD/HDE) 

tU, 

(A) 

<552, 

(A) 

<U, 

(A) 

6A2, 

(A) 

d* 

X 

(A) 

iP, 

(A) 

W 

(A) 

iW 

(A) 

12323 

15 

0.5 

0.5 

05 

068 

1548.0 

05 

0.7 

15 

08 

13621 

4.0 

0.2 

0.4 

0.4 

076 

1547.7 

27 

05 

3.8 

08 

21856 

1J 

0.5 

1.0 

05 

065 

15507 

02 

0.3 

15 

08 

24131 

12 

0-5 

0.9 

0.5 

069 

1548.0 

02 

03 

1.7 

08 

34816 

3.1 

0.8 

0.4 

05 

079 

1547.6 

21 

07 

3.4 

08 

35299 

0.1 

0.4 

0.4 

04 

074 

1548-2 

-02 

06 

08 

05 

35653 

7.0 

1.0 

as 

04 

082 

1545.4 

5.1 

1.0 

62 

08 

36430* 

... 

... 

... 


... 

... 

ao 

... 

00 


36822 

6.4 

1.0 

12 

05 

0.75 

1547.7 

3.9 

1.0 

4.7 

08 

36824 

2.9 

0.4 

0.8 

04 

037 

1550.9 

05 

0.8 

08 

08 

37209 

0.9 

0-3 

0.8 

03 

065 

1548.2 

0.1 

08 

1.4 

0.8 

37366 

65 

0.4 

0.5 

04 

070 

15482 

45 

0.6 

4.4 

08 

66665 

4.1 

0.4 

1.6 

04 

098 

1548.1 

25 

08 

5.8 

0.8 

120307 

0.8 

0.4 

0.9 

04 

0.52 

1550.8 

-0.1 

03 

0.9 

08 

132058 

1.1 

0.4 

0.6 

04 

0.65 

1547.9 

03 

04 

12 

08 

142669 

1.1 

0-3 

0.4 

03 

039 

1547.8 

05 

05 

07 

0.8 

143275 

3.3 

1.0 

05 

04 

065 

15475 

1.9 

0.9 

29 

08 

144470 

0.9 

0-5 

0.8 

05 

050 

1550.7 

01 

05 

15 

0.8 

148703 

1.1 

0.4 

0.6 

0.4 

053 

15482 

05 

0.5 

1.1 

0.8 

149438 

6.5 

0.6 

1.7 

0.4 

0.85 

1547.5 

4.1 

0.8 

5.8 

0.8 

151985 

1.0 

0.4 

0.9 

0.4 

051 

15482 

ai 

05 

u 

0.8 

152246 

13.9 

1.0 

0.9 

05 

0.87 

15421 

115 

15 

10.7 

0.8 

163892 

8.5 

0.8 

0.7 

0.5 

0.87 

1540.9 

6.8 

1.1 

7.9 

0.8 

201345 

29 

0.5 

tt7 

05 

078 

1550.0 

1.7 

0.7 

3.7 

0.8 

207538 

6.6 

0.5 

0.6 

05 

0.77 

1546.4 

4.6 

0.7 

6.2 

0.8 

208440 

12 

0.5 

0.6 

05 

0.63 

1550.4 

0.4 

0.6 

1.6 

0.8 

214680 

6.8 

0.5 

0.9 

05 

078 

1546.3 

4.7 

08 

62 

0.8 

217101 

1-3 

0l5 

0.8 

05 

045 

15505 

02 

0.4 

1.0 

0.8 

218195 

127 

12 

-0-2 

08 

0.96 

1536.9 

124 

1.9 

10.0 

08 

218342 

7.0 

0.5 

08 

05 

073 

15475 

45 

0.7 

5.4 

08 

251847 

3 2 

0.5 

08 

05 

081 

1548.0 

1.9 

0.7 

3.7 

08 


* Uncertainty in d: 0.05 in continuum units. 

* i 4 = wavelength of maximum profile depth. 

* Line probably not present. 



Fig. 8. — Dependence of the wind parameter on effective temperature. Stars that are labeled are discussed further in the test The symbols are the same as in Fig. 1. 
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(Paper I). It is dear from Figure 8 that the temperature where 
asymmetries become measurabJe for stars of lower gravity is 
between 22,000 and 25,000 K, but with our limited data on 
such stars we are unable to constrain “ tum-on " temperatures 
further. 

Another interesting aspect of Figure 8 is the position of the 
“enhanced wind” star (Walbom et aL 1985) r Sco (HD 
149438). Its position in that figure does not imply that its wind 
is enhanced for its effective temperature; on the contrary, com- 
pared to stars with effective temperatures near 30,000 K, its P w 
is somewhat smaller. However, t Sco does possess a signifi- 
cantly larger surface gravity than those stars. If it were com- 


pared to stars of its gravity, its wind might still turn out to be 
enhanced. 

In spite of having nearly the same value of P w as the three 
stars at 30,000 K, r Sco has a C rv profile that is markedly 
different: much less asymmetric and less shifted to shorter 
wavelengths, but significantly deeper (see Fig. 9a). In view of 
the strength of this star's UV resonance lines, it is interesting 
that Waters et aL (1993) have tentatively identified this star as a 
pole-on Be star based upon emission in certain hydrogen 
Brackett lines. This connection should be investigated farther. 

Like t Sco, HD 66665 has a nearly normal value of P w but 
unusually symmetrical, deep C rv resonance lines. In fact, these 



Wavelength (A) 

Fig. 9 a 



Wavelength (A) 

Fig. 96 

Fig. 9.— Comparison of the C rv profiles in selected stars. Note the similarity in the profiles of t Sco (T*, =■ 31,500 It, log g — 4.25) and HD 66665 (7^, — 28,000 
K, log g — 4.0). The profile of 10 Lac (7^, =* 30,000 K, log g =» 4.0) is markedly different from that of r Sco. 
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lines greatly resemble those of r Sco, with the exception that 
they appear even more photospheric and deeper (Figs. 9 a and 
9b). This star has a projected rotational velocity that is among 
the smallest in the sample and possesses by far the strongest 
photospheric C iv resonance lines of any star in the program; 
as previously noted, other UV resonance lines are also 
unusually strong (Grigsby 1992). Because of this star’s simi- 
larity to r Sco, observations of its H Brackett lines would be of 
great interest. 

Figure 10 shows the temperature dependence of C rv equiva- 
lent width as measured by the hand-drawn profiles. Superim- 
posed upon this figure are synthetic photospheric C tv 
equivalent widths produced by PAM (Paper I). The similarity 
between the trends of P w and equivalent width with tem- 
perature suggests that P w is related to overall carbon absorp- 
tion. The one point that is clearly different is that representing 
HD 66665. The great strength of the C rv doublet in this star is 
clearly shown in Figure 10. 

The dependences of the wind parameter and C rv equivalent 
width on luminosity are shown in Figures 11a and lib; we 
offer these plots as observational results. Both P w and equiva- 
lent width show a strong correlation with luminosity. An 
attempted eyeball fit of a straight line to the data in Figure 1 la 
suggests that there is a sudden drop in P„ at a luminosity 
somewhere between 4.2 and 4.5 dex. Clearly, the C rv doublet 
does not show evidence of winds for log (L/L 0 ) < 4.0. 

The apparent horizontal clustering of the data points may 
not be significant, given the size of the error bars: the 
maximum scatter is about 2 a. It has been suggested (see the 
review by Bedding 1993) that & Sco (HD 143275) has an anom- 
alously low mass-loss rate for its luminosity. In Figure 11a, 
however, it is not at all unusual. It does lie below the main 
trend in Figure lib, but HDE 251847 lies very close to it in 
that figure. For both stars, the lower extension of the error bar 
in luminosity would place them in the zone where P„ drops 
steeply to zero. In short, both are probably transition objects. 
Much more clearly anomlalous for its luminosity is HD 13621. 


This star, is, however, normal for its temperature and gravity. 
In view of its unusually high luminosity for its temperature and 
gravity (Fig 1), perhaps the luminosity is the source of the 
anomaly. The possibility that this star may have an undis- 
covered, equally luminous companion seems unlikely because 
our model atmospheres match the optical spectrum unusually 
well (Paper I)- Perhaps HD 13621 lies closer to us than the 
association. 

Our sample has seven stars in common with those of 
Howarth & Prinja (1989) and Prinja (1989), who used line- 
profile fitting to determine the product of the mass-loss rate 
and the ionization fraction of C"' 3 . Figure 12 shows a tight 
correlation between this quantity and P„. This result supports 
our conclusion that P w is a good stellar wind indicator for 
chemically normal stars. Both our Figure 11a and Figure 8 
from Prinja (1989) are consistent with the hypothesis that mass 
loss “turns on” suddenly in main-sequence stars at log 
(L/L 0 ) ~ 4.5 and strengthens steeply with increasing lumi- 
nosity above that value. Of course, the steep slope just above 
the threshold could be caused by a rapid increase in the ioniza- 
tion fraction of C +3 , and the slope of log M itself could be a 
constant. 

6. COMPARISON WITH STATIC, NON-LTE MODEL 
PHOTOSPHERES 

Figure 8 shows that, in near-main-sequence stars, the C rv 
absorption lines acquire their asymmetry and cease to be 
purely photospheric in origin at effective temperatures between 
22,000 and 27,000 K. Although we did not compute model 
stellar winds, we attempted to gain insight into the cause of 
asymmetrical C rv absorption by studying the physical condi- 
tions related to C rv line formation in our model photospheres. 
The good agreement between the synthetic and the observed 
C rv equivalent widths at T M =» 25,000 K (§ 5 and Fig 10) 
suggests that PAM simulates these conditions realistically. 

The temperature structure of two models is shown in Figure 
13, and the figure also illustrates the relationship between 



Fio. 10. — Temperature dependence of C rv equivalent width in the program stars. Diamonds represent mode! photospheric equivalent widths as predicted by the 
tioos of Paper 1. Other symbols are as in Fig. 1. 
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Fig. 116 


Fig. 1 1. — Dependence of (a) wind parameter and (6) C rv equivalent width as a function of luminosity (in solar units) 


column density and level number. The continuum emerges 
from roughly level 40, a region several thousand degrees hotter 
than the line-forming region. Figure 14 shows the PAM predic- 
tions of the populations of the ground state of C rv (relative to 
all other carbon atomic states; the sum of the populations of 
all carbon states is unity) as a function of effective temperature 
for five levels in the model atmosphere. Strikingly, at effective 
temperatures above 27,000 K, the C rv ground-state popu- 
lation is larger at level 1 than at level 40. This is a completely 
non-LTE result ! At 31,500 K the population of level 10 also 
exceeds that of level 40. The populations at levels 20 and 30 are 
depressed relative to those at level 40 in all models; these 
populations are typically smaller than those at level 1 by an 


order of magnitude or more. The model also indicates that the 
line core optical depth of the shortward component has not 
begun to approach unity in level 1 at an effective temperature 
of 25,000 K, but it is roughly unity in level 1 at 28,000 K, and at 
an effective temperature of 30,000 K it is 16. 

7. DISCUSSION 

These model predictions thus indicate that, at about 
—28,000 K, the C rv line cores become optically thick at the 
top of the photosphere while remaining optically thin in the 
intermediate layers between the top and the level of continuum 
formation. Therefore, it is tempting to hypothesize that radi- 
ation pressure on C rv actually begins to drive the wind at this 



808 


GRIGSBY & MORRISON 


Vol. 442 



Wind Parameter (A) 

Fig. 12. — Correlation of the product of the mass-loss rate and C*’ ionization fraction (Howarth & Prinja 1989rPrinja 1989) with P„, in the seven stars they had 
in common with the present work. 


temperature. Both the opacity enhancement at the top of the 
photosphere and the relative transparency below should con- 
tribute to radiation pressure on C rv at the top of the photo- 
sphere, which could lead to the establishment of a wind. 

Of course, radiation pressure on C rv ions does not explain 
cases like HD 12323, which shows a wind signature in the N v 
lines but not in C rv. In such C-poor stars, perhaps the N v 
lines play a similar role in initiating observable mass loss at a 
similar temperature. A full explanation of this effect should 
include line transfer in a moving atmosphere, preferably by 
means of a unified model atmosphere. 

The alternative, more conventional, explanation for the 
transition from symmetrical to asymmetric C rv line profiles is 


that radiation-driven winds exist at temperatures and lumi- 
nosities below this transition but are unobservable because the 
optical depth in the wind in the C rv doublet is too small. The 
transition from unobservability to observability would then be 
caused by an increase in the line core opacity of the outermost, 
moving layers of the atmosphere. The non-LTE population 
increase at the top of the photosphere, described above, could 
also explain such an opacity effect if it continued outward into 
the wind. 

Therefore, our model photospheres are not capable of dis- 
criminating between these two hypotheses. In the following 
paragraphs, we will consider other arguments. 

It is not implausible that the C tv resonance doublet could 



Log (m) g cm' 2 

Fig. 13. — Electron temperature as a function of column density in two PAM models. The solid vertical lines illustrate the relationship between model level and 
column density. 
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effective Temperature 

FiO. 14. — Model prediction for the relative population of the C rv ground state as a function of effective temperature. The point representing level 1 at 
■» 23,000 K was unavailable and is missing. The ground-state population at level 1 (top of photosphere) exceeds that of level 40 at temperatures higher than 
27,000 ft. 

be an important wind driver since it is the strongest unsatu- correction. According to Friend & Abbott (1986), the effect of 

rated line near the peak of the energy distribution at this effec- that correction is to decrease the line acceleration at the base of 

tive temperature. In the temperature regime under study, C, N, the wind by a factor that can be as large as 0.6. With a decrease 

and O contribute more of the line acceleration than any other by this factor, we estimate that the static limit would move up 

group of elements (Abbott 1982). When the wind is weak (small in the H-R diagram by about A log L- 0.2. Although this 

optical depth), the strongest lines are barely desaturated, and change would improve agreement with our observed limit, 

those lines dominate the radiation force. This situation, briefly caution is in order because the magnitude of this shift is corn- 

mentioned in Abbott (1982), applies here. parable to observational error and because other approx- 

Another argument in favor of the hypothesis that radiation- imations in Abbott’s theory still need to be removed. For this 

driven mass loss occurs only where it is actually observed in reason, a detailed comparison between the static limit as calcu- 

the C rv asymmetry is the general agreement between the loca- lated by Abbott (1979) and the boundary of the domain of 

tion on the main sequence where winds appear and the loca- observability of C rv profile asymmetry found in this study 

tion in the H-R diagram of the “ static limit " computed by would be premature. It would be of great interest for the theo- 

Abbott (1979) on the basis of radiation-driven wind theory. retical static limit to be recalculated by more rigorous methods 

The static limit is the locus of points in the H-R diagram at (Kudritzki et al. 1989). 

which the line acceleration is equal to the effective acceleration Perhaps the most sensitive diagnostic of a radiation-driven 

of gravity. Above this locus, static atmospheres cannot exist wind is emission of X-rays, which are believed to originate in 
and radiation-driven winds are inevitable. Below the locus, shocks caused by instabilities in the wind. According to cur- 
radiation pressure can sustain a wind if another mechanism rently available evidence. X-ray emission persists to lower T e(r 

sets up a velocity gradient in the atmosphere. Therefore, one on the main sequence than does C rv profile asymmetry. For 

expects the static limit to correspond, not necessarily to a example, Cassinelli et al. (1994) found that main-sequence stars 

boundary in the H-R diagram where winds disappear com- as low in effective temperature as 14,000 K. show definite X-ray 
pletely, but at least to a boundary where the incidence of stellar emission. On its face, this result is strong evidence that the 

wind drops. Various authors, such as Abbott (1985b), have threshold in C rv profile asymmetry is not a real wind thresh- 

pointed out the general agreement between the theoretical old but only an optical depth effect. 

static limit and the observational lower luminosity limit for But a confusing factor in the statistics of X-ray emission in 

observation of winds in ultraviolet resonance lines in normal hot stars is that binaries with companions of unknown spectral 
stars. type are present in the sample. Coronal X-ray emission from a 

On the other hand, the relevance of the static limit can be late-type companion would account for the observed flux of 

questioned. Perhaps many or all early B-type stars, by means X-rays from mid-B-type stars in binaries (Meurs et al. 1992). In 

of rapid rotation or pulsation, set up enough of a velocity the sample of Cassinelli et al. (1994), only one of the three 

gradient for radiation pressure to sustain a wind. In addition, coolest stars detected (14,400-18,000 K) is not a binary with a 

the precise location of the static limit is not well determined by companion of unknown type. Still, that sample, though admit- 

existing theory. Abbott’s (1979) pioneering calcuation made tedly small and heterogeneous, includes several apparently 

several approximations, one of which is that the star is a point single. X-ray-emitting stars that are cooler than 26,000 K, and 

source as seen from the wind. More recent calculations have it is difficult to explain all of them away, 

removed this approximation by incorporating the finite disk X-rays in a larger sample of early B-type stars were studied 
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by Meurs et al. (1992). In their Figure 3 is presented a histo- 
gram of the fraction of stars in which X-rays were detected with 
ROSAT as a function of spectral type. The fraction decreases 
strongly toward later type until it reaches a background level 
at which, on average, — 10% of the stars in each spectral class 
were detected. Although other explanations are possible, it 
may be that these background detections are binary systems 
that include coronally emitting, late-type components. What- 
ever its origin, this background is probably not caused by 
normal, radiation-driven winds. The latest spectral type at 
which the X-ray detection rate is more than 1 a above this 
background is B1 ; a 2 a criterion would place the threshold at 
BO. Therefore, it is not yet possible to rule out the hypothesis 
that the detectability threshold for radiation-driven winds in 
hot stars is at the same location in the H-R diagram for X-rays 
as it is for profile asymmetries in the C rv resonance doublet 

MacFarlane, Cohen, & Wang (1994) have shown that X-ray 
emission at the observed levels, if it originates in the stellar 
wind, significantly affects the ionization fractions in the winds 
and perhaps even in the photospheres of hot stars. At 
25,000 K, for example, they find an enhancement in the ioniza- 
tion fraction of C rv by at least four orders of magnitude in the 
wind (at v = 0.50 compared to a model without X-rays, and 
they suggest that the X-ray emission results in a major 
enhancement of highly ionized species at the base of the wind 
also. In our non-LTE model photospheres, however, the com- 
puted and the observed C rv equivalent widths agree very well 
at and below an effective temperature of 28,000 K, without 
inclusion of X-rays or any other “superionization" mecha- 
nism. This result implies that X-rays are unimportant in the 
ionization balance of C rv in the photospheres of early B-type 
stars. 

8. CONCLUSIONS 

From our examination of the morphology of the C rv 
211548, 1550 resonance doublet, we conclude that the doublet 
is first detectably asymmetric in near-main-sequence stars 
(those with log g s 4 ± 0.1) of normal abundances at effective 
temperatures of 26,000-28,000 K and log ( L/L 0 ) s 4.5, when 
the model photospheric equivalent width of the C rv doublet is 
24-3.0 A. We find little evidence that iron or other weak 
photospheric blends mimic mass loss in marginal cases; the 
most straightforward explanation for the asymmetries is that 
they are caused by weak winds. The winds make their appear- 
ance at lower effective temperatures for stars of-loWer surface 
gravity, though with our limited data we cannot identify the 
threshold effective temperature in those stars. Some evidence 
for winds is seen in the N v and Si rv resonance doublets in 
some of the program stars of normal composition that possess 
evidence of winds, but the mass-loss signature is much stronger 
in the C rv doublet. We see little evidence for the existence of 
“enhanced mass loss " stars as described by other authors (e.g., 
Walbom et al. 1985), but stars such as r Sco and HD 66665 do 
present unique C rv profiles and demand additional study. The 
problem of similarly classified stars with different C iv mor- 


phology that has been found by other workers appears par- 
tially solved, presumably because of the quantitative 
measurement of effective temperature and surface gravity in 
our program stars. Nevertheless, we still see cases in which 
stars with similar effective temperatures and surface gravities 
display markedly different C iv profiles (e.g^ HD 13621, HD 
34816, HD 66665, and HDE 251847 at an effective temperature 
of 28,000 K). 

The initial appearance of the wind in this doublet coincides 
in effective temperature with a non-LTE effect: the population 
of the model C tv ground state is larger at the top of the 
photosphere than it is at the level of continuum formation for 
effective temperatures higher than 27,000 K; at intermediate 
photospheric levels, the population of the C rv ground state is 
an order of magnitude or more smaller than it is at the top of 
the photosphere. Non-LTE models also predict that the short- 
ward component of the C rv doublet reaches optical depth 
unity at a column density of —5 dex. This combination of 
effects allows the underlying photosphere to exert a maximum 
radiative force on the top layers and may initiate mass motion 
in the photosphere. At 28,000 K the asymmetry in the doublet 
is measurable, and at higher temperatures the existence of a 
wind is evident through its asymmetry and blueshift. This sce- 
nario does not account for winds seen in the N v doublet in 
N-rich stars, which are carbon poor; that problem needs 
further investigation. 

Section 7 of the paper concluded with a discussion of 
whether the boundary we have identified in the H-R diagram 
represented a real onset of radiation-driven winds or merely a 
threshold in the optical depth of the winds in the C rv reso- 
nance doublet. Neither the location of the theoretical “static 
limit " nor the incidence of X-ray emission in early B-type stars 
can yet be used to rule out the possibility that the threshold is 
real. Further work, both observational and theoretical, is 
needed to resolve this question. 
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